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ABSTRACT 



Context. The detection of pulsational frequencies in stellar photometry is required as input for asteroseismological modelling. The 

second short ran (SRa02) of the CoRoT mission has provided photometric data of unprecedented quality and time-coverage for a 

number of O-type stars. 

Aims. We analyse the CoRoT data corresponding to three hot O-type stars, describing the properties of their light curves and we 

search for pulsational frequencies, which we then compare to theoretical model predictions. 

Methods. We determine the amplitude spectrum of the data, using the Lomb-Scargle and a multifrequency HMM-like technique. 

Frequencies are extracted by prewhitening, and their significance is evaluated under the assumption that the light curve is dominated 

by red noise. We search for harmonics, linear combinations and regular spacings among these frequencies. We use simulations with 

the same time sampling as the data as a powerful tool to judge the significance of our results. From the theoretical point of view, we 

use the MAD non-adiabatic pulsation code to determine the expected frequencies of excited modes. 

Results. A substantial number of frequencies is listed, but none can be convincingly identified as being connected to pulsations. The 

amplitude spectrum is dominated by red noise. Theoretical modelling shows that all three O-type stars can have excited modes but 

the relation between the theoretical frequencies and the observed spectrum is not obvious. 

Conclusions. The dominant red noise component in the hot O-type stars studied here clearly points to a different origin than the 

pulsations seen in cooler O stars. The physical cause of this red noise is unclear, but we speculate on the possibility of sub-surface 

convection, granulation, or stellar wind inhomogeneities being responsible. 

Key words, stars: variables: general - stars: early-type - stars: individual: HD 46223 - stars: individual: HD 46150 - stars: individual: 
HD 46966 - stars: oscillations 



1. Introduction 

Even though the O-type stars are located in the Hertzsprung- 
Russell (HR) diagram inside a zone where pulsations are ex- 
pected, no asteroseismological modelling of these stars could be 
performed prior to the CoRoT (Convection, Rotation and plane- 
tary Transits, Bagli n et al.|2006 Auvergne et al. 2009) mission. 
Actually, few examples of variability have been identified. The 
most obvious ones are the 09.5 V stars C, Oph (Kambe et al. 
[1997) and HD 93521 ( [Howarth & Reid|1993[|Rauw et al.|2008E 
and their variability is likely related to non-radial pulsations with 
periods of a few hours. However, the majority of these detections 
were made spectroscopically. The reason is that the amplitudes 
of pulsations are too low to be measured from ground-based pho- 



tometry. Furthermore, these pulsations can be contaminated by 
variable stellar winds. 

The second short run (SRa02) of the CoRoT satellite in the 
Asteroseismologic Channel was partly devoted to the investiga- 
tion of the photometric variability of O-type stars. Pointing to- 
wards the anti-centre of the Galaxy, this instrument observed ob- 
jects belonging to the young open cluster NGC 2244 inside the 
Rosette nebula and to the surrounding association Mon OB2. 
The CoRoT data of three O-type stars that are part of this run 



have been analysed in previous papers: HD 46149 by Degroote 



etaqpOlO) HP 47 129 by |Mahy et al.lpOTT) and HD 46202 by 



Briquet et al. (201 1 ), the latter paper containing the only forward 



* The CoRoT space mission was developed and is operated by the 
French space agency CNES, with participation of ESA's RSSD and 
Science Programmes, Austria, Belgium, Brazil, Germany and Spain. 
** Tables 2, 3 and 4 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79. 128. 5) or 



modelling in terms of seismic data of an O star so far. 

In the present paper, we study the three remaining O-type 
stars that are part of the SRa02 run. One of the targets is 
HD 46223, t he hottest member of NG C 2244. Situated at about 
1.4—1.7 kpc (Hensberge et al. 20001, this star was previously 
quoted in the literature as having a spectral type 05 V ( Hiltner 
[1956] IBisiacchi et al.||T982), 05 ((f)) (|Conti & Leep ^ 



via http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 
*** Senior Research Associate FR.S.-FNRS 



04 ((f)) (Masseyetal.[1995) and 04 V((f + )) (Mahyetal. 



19747 



2009) 



The 9-year spectroscopic campaign by |Mahy et al.| to investi- 
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2. HD 46223 

2.1. CoRoTdata 

The data we use here were obtained during the CoRoT second 
short run (SRa02) mad e from 08 Oct to 12 Nov 2008. We start 
from the level N2 data (Sam adi et al.|20 07), which for each ob- 
servation list the time, flux, error on the flux and flagging in- 
formation. The flagging indicates whether these data are con- 
taminated by instrumental and environmental conditions of the 
CoRoT satellite such as, e.g., the South Atlantic Anomaly (SAA) 
and other Earth orbit perturbations (see Auvergne et al. 2009 1. In 



our analysis, we discard all flagged points, as well as those few 
points that are listed with a negative flux error. 



ILx.jaJjjJ 2.2. Spectral window 
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Fig. 1. Spectral window of the HD 46223 observations. The inset 
zooms in on the low-frequency region. 



gate multiplicity did not reveal any significant variation of the ra- 
dial velocity related to the presence of a companion. In addition, 
Chandra mosaic observations of NGC 2244 ( [Wang et aL|2 008 ) 
detect no strong clustering of X-ray sources around HD 46223, 
leading to the conclusion that there is also no lower-mass neigh- 
bour present. According to Wang et al. the small number of 



close-by X-ray sources could indicate that this star is younger 
in comparison to the population of the central part of the cluster. 
The next target is HD 46150, the second hottest star in 
NGC 2244. Literature values for the spectral type are between 
05V((f)) ( |Underhill & Gilroy| [1990} IMaiz-Apellaniz et al. 
2004] and Q5.5 V((f)) ( |Conti & Ebbets||19771 |Garmany et al. 
1980 1. The detailed study of |Mahy et al.| ( |2009| l confirms the 
05.5 V((f)) spectral type. Radial velocity variations have been 
detected in this star, but it is not clear if these are due to binarity 
or motions in the stellar atmosphere (see Mahy et al. 2009 , and 
references therein). Mahy et al. monitored this star but could not 



find any significant period in their data. They note the lack of 
spectral variability over short timescales (~ 7 days). 

The third target studied in the present paper is the late O- 
type star HD 46966, belonging to the Mon OB2 association. The 
spectral classification gen erally a grees on Q8 .5 V ( |Conti et aL 
T977)|Garmany et al.|198 ; Mahy et al.|2009|l or Q8 V ( jMunan 
& Tomasella 1999 Maiz-Apellaniz et al. 2004 1. Moreover, spec- 
troscopic monitoring reveals no clues on the existence of a sec- 
ondary component (Ma hy et al.|2009) . 

These three stars are among the youngest (1-6 Myr, B onatto| 
& Bica 2009) and most massive observed by the CoRoT satellite. 
They therefore have the potential to extend our knowledge about 
pulsations into the higher-mass domain. The CoRoT light curves 
are of unprecedented quality and have a long observational time 
span (longer than the rotational cycle). They thus provide the 
ideal data set to reveal the presence or absence of such pulsa- 
tions. 

We first present, in Sect. [2] the detailed frequency analysis of 
the HD 46223 CoRoT light curve. In Sect.|3]and Sect.|4] we give 
the frequency analysis of HD 46150 and HD 46966 respectively. 
Section B] is devoted to determining the theoretically expected 
pulsation frequencies. We discuss the results in Sect. [6] and in 
Sect. [7] we provide the conclusions of our research. 



The spectral window for the HD 46223 observations below 
10 d _1 (Fig. lb is characterised by a first peak at / = 2.007 d _1 
and its double at about / = 4.01 1 d _1 . This is related to the pas- 
sage of the satellite, twice in a sidereal day, through the SAA. 
The relative amplitude of these two peaks is 3.5 % and 1.1 %, 
respectively. In addition, the gaps in the data set, due to the or- 
bital period of the satellite (6184 s), generate other structures 
with peaks around / = 13.972 d _1 and their harmonics. 

The spectral window also exhibits a peak at / = 2699.76 d _1 , 
corresponding to the sampling regularity of about 32.003 s. As a 
consequence, the pseudo-Nyquist frequency is located at about 
1350 d . The word "pseudo" emphasizes the fact that the sam- 
pling is not ideally regular because of the flagged data. The duty 
cycle of the observations is 86.75 %. Because the data for all 
three stars were collected simultaneously with the same instru- 
ment, the spectral window properties for the other stars are very 
similar (see Table[T|. 

2.3. Defending 

In common with almost all CoRoT targets, there is a long-term 
trend visible in the data. This decreasing slope is likely due to 
the CCD ageing (Auvergne et al. 2009). To remove this trend we 
divide the flux counts by the best-fit linear slope to the data. 

The final version of the CoRoT light curve covers T ~ 
34.3 days. The exact start and end times, duration, number of 
points and frequency resolution are listed in Table [T] We do not 
convert the fluxes to magnitudes. 



2.4. Frequency-by-frequency prewhitening 

The range of variations observed in the HD 46223 CoRoT light 
curve (Fig. [2] top) is of order 16,000 counts, corresponding to 
8 mmag. No clear pattern in the variations is visible. In order 
to better understand the variability present in the CoRoT light 
curve, we apply a Fourier analysis based on the Lomb-Scargle 
method ( |Lomb|1976||Scargle|1982]l 

The amplitude spectrum (Fig. |3] top) exhibits many peaks 
and appears to be noisy. The highest peak is located at about / = 
0.75 d _1 and the second one at about / = 0.25 d , but no high- 
significance peaks are clearly detected, nor does there appear 
to be systematics, such as a constant spacing between peaks. A 
plot of the spectrum up to the pseudo-Nyquist frequency shows 
very little power at high frequencies. In the further analysis, we 
therefore limit the frequency domain to 0-100 d _1 . 

For a more detailed analysis, we perform a traditional 
prewhitening. We start by subtracting the average from the 
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Table 1. CoRoT data and results from the analysis for the three stars. 



Parameter 



HD 46223 



HD 46150 



HD 46966 



CoRoT data 

number of raw data, N R 

number of non-flagged data, N 

duty cycle (N/N R , in %) 

start date non-flagged data (HJD - 2450000) 

end date non-flagged data (HJD - 2450000) 

duration, T (d) 

freq resolution (d _1 ) 



Spectral window 

window peak (%) at 2.0 d" 1 

window peak (%) at 4.0 d~' 
Stopping criterion: number of frequencies 

AIQ. 

BIC 

HQC 



92687 


92695 


92688 


80403 


80282 


82273 


86.75 


86.61 


88.76 


4748.488547 


4748.485591 


4748.488454 


4782.819650 


4782.819657 


4782.805898 


34.331103 


34.334066 


34.317444 


0.029128 


0.029126 


0.029140 


3.5 


3.8 


3.9 


1.1 


1.0 


1.1 


860 


821 


595 


491 


444 


277 


662 


625 


427 



Number of red-noise significant frequencies 

59 


50 


53 


Number of red-noise significant frequencies present in both halves of the observing 

10 7 


run 

3 


Red noise - fit parameters 

a a (counts) 442 ± 3 1 
t (d) 0.09 ± 0.01 
y 0.96 ± 0.02 


430 ± 65 
0.08 ± 0.02 
0.96 ± 0.03 


416 ±85 
0.17 ±0.05 
0.91 ± 0.02 



fluxes. We then determine the amplitude spectrum, using the 
Lomb-Scargle method. We select the frequency corresponding 
to the highest peak in the amplitude spectrum and fit a sine 
function (Aj sin(27r/)?,' + <pj)) to the fluxes, as a function of time 
(tt). In the fitting procedure, the semi-amplitude (Ay) and phase 
{(pj) are free parameters, while the frequency (/)) is allowed to 
vary within a range of +0.1/T from its peak valuaH(where l/T 
is the frequency resolution and 0.1 /T the step size in our fre- 
quency grid). The fitted sine function is then subtracted from 
the data and the amplitude spectrum is re-computed from these 
prewhitened data. This procedure is continued until we reach the 
noise level in the observations. After the prewhitening, the ob- 
served variability, whatever its nature, is thus described (down 
to the noise level) as a sum of terms, with each term being a sine 
function at a particular frequency. 



2.5. Stopping criteria 

To judge when we have reached the noise level, we explore 
three stopping criteria that have been discussed in the literature 
( |Hannan & Quinn|1979||Hannan|1980||Liddre|20071 l. The infor- 
mation content of the fit is given by: 



c N (m) = -2ln(L N (m)) 



+ < 



2m + 



2m(m + 1) 



-.- r- Akaike, corrected (AIC C ) 

(N - m - 1) v c/ 

m ln(N) Bayesian (BIC) > W 

2m ln(ln(/V)) Hannan - Quinn (HQC) 



where N is the number of observations, m is the number of pa- 
rameters and Lfj{m) is the likelihood function. To within a con- 



1 The ±0.1/T value turns out to be quite close to (slightly higher 
than) the formal errors we get on the frequencies (see Eq. |3j|5] and 
Tables|2]3and4). 



stant, 
15.1): 



-21n(L;v(m)) can be written as ( |Press et al.||1992| Sect. 



l 

;=1 



(yi-Fi^&J) 2 



(2) 



where y, are the fluxes observed at times f,, and 0„, is the vector 
containing the m fit-parameters. For the value of <r, we take the 
error bars listed with each observed flux. We stop the prewhiten- 
ing process when caK'«) has reached a minimum. This is very 
similar to least-squares minimization, but the additional terms 
in Eq. [T] correct for the fact that a model fit improves with an 
increasing number of parameters. We set m - 3k, counting 
3 parameters (frequency, amplitude, phase) for each of the k 
sine functions we fit. Note that we subtract the average from 
the observed fluxes before starting the prewhitening procedure 
(Sect. |2~4| i, which could lead to a bias. 

Applying this to the HD 46223 data, we find that the mini- 
mum of AIC C is reached at 860 terms, BIC at 491 and HQC at 
662. To judge the reliability of these stopping criteria, we run a 
number of simulations where we construct artificial light curves 
with a given number of sine functions. We use the same time 
sampling as the observations and attribute to it the same noise 
level and flagging (this approach is used in all simulated light 
curves discussed in this paper). For the amplitudes and frequen- 
cies we take either values from the analysis of the observations, 
or randomly drawn values distributed according to the statistics 
of the observed spectrum. Noise is added to the simulation, using 
a Gaussian random variable (with mean zero and standard devi- 
ation unity) multiplied by the listed error on the observed flux. 
We then apply our prewhitening procedure to the simulations 
and compare the results of the stopping criteria with the known 
number of terms we introduced. It turns out that BIC gives the 
best value, which is always within a factor 1.5 of the correct re- 
sult. Degroote et al. (2009b) also used BIC in preference over 
the (uncorrected) AIC as a stopping criterion in their study of 
the ft Cep star HD 1 80642. 
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Fig. 2. Detrended CoRoT light curves, covering the full 34.3 days of the observing run. Each observation is represented by a dot. 
Flagged data are not plotted. Top: HD 46223, middle: HD 46150, bottom: HD 46966. 



As a further test, we also redo the analysis on the observed 
data after rebinning over a number of points (up to 32 points, i.e. 
17min), but find approximately the same results. This is consis- 
tent with the fact that there is little power at higher frequencies. 
Based on the above, we take the BIC rounded-off value of 500 
terms as the stopping criterion. The final 500 frequencies and 
corresponding semi-amplitudes and phases are listed in Table [2] 
(the full Table [2] is only available in electronic form at CDS). 
The phases are defined with respect to t = corresponding to 
the first non-flagged data point (as listed in TablefTli. 

New frequencies found during the prewhitening procedure 
can be quite close to previously found frequencies. The reso- 
lution criterion for this is that frequencies should be separated 
by at least 1.5/T to be considered unique (Loumos & Deeming 
|1978) l. Applying this to our list of 500 frequencies, we find that 
276 of them are unique. 

The canonical errors estimated on these frequencies (/), 
and their corresponding semi-am plitudes (A/) and phases (0/) 
are given by the expressions of |Lucy & Sweeney] ( |1971[ ) and 



Montgomery & Q'Donoghue| ( [T999] l: 

e(f) = 




(3) 

(4) 
(5) 



with ctj the standard deviation on the light curve at the cur- 
rent stage of prewhitening. These formal errors are also listed 
in Table |2] 



2.6. Multifrequency fitting 

The Lomb-Scargle based iterative prewhitening consists in a 
one-by-one determination of the sine functions and in their re- 
moval from the light curve. In a gapped or unevenly sampled 
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Fig. 3. Amplitude spectrum of the CoRoT light curves computed by the Lomb-Scargle method. The insets show a zoom-in of the 
amplitude spectrum on the low-frequencies domain. Top: HD 46223, middle: HD 46150, bottom: HD 46966. 



data set, the height of peaks is dependent on the height of other 
peaks in the periodogram. A more general procedure should 
therefore determine all sine function terms simultaneously. 

In order to be more reliable in the detection of the fre- 
quencies, we also use a multifrequency algorithm introduced by 
|Gosset et aLl([2001| th eir Eq. A1 3 to A1 9) and developed from 
the |Heck et al.| ( |1985] hereafter |HMM| l method. We note that 
this multifrequency algorithm is not based on the Lomb-Scargle 
technique. However, in the case of CoRoT data, we find that 
the amplitude spectrum computed from Lomb-Scargle and from 
HMM are in very good agreement. The multifrequency algo- 
rithm takes into account the mutual influence of peaks by fitting 
them all together, and is therefore in principle an improvement 
on the one-by-one prewhitening procedure used in Sect. |2.4| 
Similarly to the prewhitening procedure, the present method re- 
fines the value of each frequency (allowing it to vary within the 
natural width of the individual peaks). 

We apply this method to the non-flagged data of the 
HD 46223 CoRoT light curves but the great number of data 



points involves an excessive computation time. Accordingly, we 
have to deal with a limited number of terms at the same time. 
We find that the values of the frequencies fitted by the multifre- 
quency programme are very similar to those determined by the 
prewhitening procedure and only small differences in the ampli- 
tudes are found. We therefore do not list these results. 



2. 7. Significance levels 



Although the quality of the CoRoT curve is unparalleled, it is 
still important to check whether peaks detected in the ampli- 
tude spectrum might result from a random variation rather than 
representing a periodic signal. To do this, we use the statistical 
criterion presented by Gosset (2007) and quoted in Mahy et al. 
(2011 1. The probability that at least one of the semi-amplitudes 
of an extended set of frequencies exceeds a threshold z under the 
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Table 2. List of sine function terms for HD 46223. 



ID 


Freq. 


Freq. 


Semi-ampl. 


Semi-ampl. 


Phase 


Phase 


Standard 


Sig. 


Red- 


Half- 






Error 




Error 




Error 


Deviation 


Red- 


noise 


run 




(d" 1 ) 


(cT 1 ) 


(counts) 


(counts) 


(rad) 


(rad) 


(counts) 


noise 


Significant 


Test 


1 


0.7548 


0.0050 


989.095 


305.732 


-2.500 


0.309 


3296.681 


1.0000 




* 


2 


0.2493 


0.0068 


882.119 


371.151 


-1.830 


0.421 


3221.288 


1.0000 




* 


3 


0.5102 


0.0067 


772.225 


323.271 


1.476 


0.419 


3160.377 


1.0000 




* 


4 


0.6858 


0.0063 


752.212 


296.426 


-0.497 


0.394 


3112.720 


1.0000 




* 


5 


1.1944 


0.0056 


699.589 


244.181 


0.676 


0.349 


3067.088 


1.0000 




* 


6 


1.8456 


0.0046 


691.462 


199.852 


-2.537 


0.289 


3026.889 


1.0000 




* 


7 


0.0395 


0.0091 


681.390 


385.046 


-1.697 


0.565 


2986.975 


1.0000 




* 


8 


0.4251 


0.0077 


650.217 


312.906 


-0.054 


0.481 


2948.151 


1.0000 




* 


9 


0.9477 


0.0064 


636.137 


251.785 


-1.726 


0.396 


2912.468 


1.0000 




* 


10 


1.4370 


0.0053 


642.207 


212.699 


-1.822 


0.331 


2877.663 


1.0000 




* 


59 


7.5519 


0.0027 


303.826 


51.957 


0.500 


0.171 


1913.541 


0.0080 


R 





Notes. This is an excerpt of Table 2. The full table is only available in electronic form at the CDS via anonymous ftp to cdsarc . u-strasbg . fr 
(130. 79. 128. 5) or via http : //cdsweb . u-strasbg. fr/cgi-bin/qcat?J/A+A/ We list the first ten terms as well as the first term that is 
significant under the assumption of red noise (indicated with 'R'). Terms that are found in both halves of the observing run are indicated with '*' 
in the last column. 



null hypothesis of a stochastic process of variance cr 2 f is given 
by 



Prob[Z max > z] = 1 - e 



,-0.93z + ln(0.8A0 



(6) 



where Z max = maxy<y Ny Z(/), N is the number of points in the 
data set and Z(/) = A 2 f N/(4cr 2 f ) with Aj the semi-amplitude of 
the frequency / and o~f the standard deviation on the light curve 
at the current stage of prewhitening. 

We also explore a different version of this criterion, which 
uses (|Degroote et al. 2009a): 

Prob[Z max > Z ] = 1 - (1 - exp(-z)f-, (7) 

where Nt is the number of independent frequencies. De grootel 
et al. show that the CoRoT data are sufficiently well equally- 

2N. In our work, the first criterion 



spaced that we can take N, = 2N. In our work. 

(given by Eq. [6]) turns out to be slightly stricter, and we therefore 

discuss only results obtained with that one. 

In applying this criterion, the simplest approach is to assume 
white noise, i.e. crj is independent of frequency and its value 
is given by the standard deviation on the light curve at the cur- 
rent stage of prewhitening. Under that assumption all 500 sine 
function terms are highly significant. 

2.8. Red noise 

It is quite clear from the amplitude spectrum (Fig. [3} that the 
power increases substantially towards lower frequencies. This 
strongly suggests the presence of red noise. We stress that the 
word "noise" does not imply an instrumental origin but rather 
points to predominantly stochastic behaviour: in Sect, 
discuss the possible physical cause of the red noise. 

In the present section, we determine the observed properties 

( 2002) 1 suggest to describe the 



6.2 



we 



of this red noise. Stanishev et al. 



amplitude spectrum a(f) by using the following fit-function: 



a{f) = 



a {) 



1 + (2mf)y 



(8) 



with ofo the scaling factor, y the slope of the linear part (in a 
log-log plot) and t an estimation of the mean duration of the 
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Fig. 4. Fit-function (Eq.^l red line) to the amplitude spectra as a 
function of frequency, in a log-log plot. Top: HD 46223, middle: 
HD 46150, bottom: HD 46966. 
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In the prewhitening procedure (Sect. |2.4| i we then need to fit a 
sine function with this frequency to the light curve, and subtract 
it. The fitting procedure will not succeed, because in attempting 
to fit only a minor contribution to the light curve, it will be led 
astray by the much higher-amplitude contributions (which have 
not been removed yet). It will therefore give an incorrect result, 
or will fail to converge. 

Applying this procedure, we find 59 significant sine func- 
tions with a significance level less than 0.01; we indicate them 
with an 'R' in Table |2] 



2.9. Split in two halves 

We next apply an additional method to determine the reliability 
of the sine function terms we found. We split the data set into 
two halves, and redo the analysis on each half separately. We 
then check how many of the original terms can also be found 
in both halves. We consider two frequencies (f\ and /i) to be 
the same if \f\ - /2I < -\Je 2 (fi) + e 2 (/2), where the errors are 
defined in Eq. [3] Simulations with artificial data confirm that 
this is the appropriate criterion to use. We could of course relax 
it by considering equality to within 2 or 3 times the 1-sigma 
error bar, but this would considerably increase the number of 
false coincidences. 

If we first look at the terms determined assuming white 
noise, we find that, of the 500 terms extracted from the full run, 
202 are present in both halves. This is substantially less than 
the 282 terms we find if we construct an artificial, noisy, light 
curve from the first 100 observed terms of HD 46223. Other 
simulations with frequencies drawn randomly according to the 
observed distribution also find a significantly larger number of 
terms in both halves. 

If we limit ourselves to the 59 red-noise significant terms 
listed in Table [2] we find that only 10 are present in both halves 
of the data set. 



Fig. 5. Autocorrelation functions (black curve) of the power peri- 
odograms and number of components in a chain with a given fre- 
quency spacing (red curve). Top: HD 46223, middle: HD 46150, 
bottom: HD 46966. 



dominant structures in the light curve. We fit this profile to the 
amplitude spectrum of HD 46223 (Fig. HI top) in order to quan- 
tify the information about the red noise contained in the peri- 
odogram. We apply a least-squares fit to the logarithm of both 
frequency and semi-amplitude, limiting the frequencies used to 
0-100 d . The resulting fit parameters and their 1-sigma error 
bars are listed in Table Q] 

To re-evaluate the significance of the peaks under the hypoth- 
esis of red noise, we need to determine its <x/; re d- We do so by 
first converting the fit-function for amplitude (Eq. [8]) into power. 
We then re-scale that function, so that its integral over frequency 
(up to the pseudo-Nyquist frequency of ~ 1350 d _1 ) is the same 
as the variance on the flux multiplied by the pseudo-Nyquist fre- 
quency. 

We then redo the prewhitening procedure and for each fre- 
quency we determine the significance of the peak by using Eq.[6] 
with Z(f) - A 2 f N/(4o- 2 fKd ). It is important to realise that in doing 
so, we still select the highest peak at each step in the prewhiten- 
ing process, not the one that is most significant under red-noise 
conditions. If we selected the most significant peak, it would 
correspond to a small-amplitude contribution to the light curve. 



2.10. Frequency combinations and spacings 

We next search for the presence of low-order harmonics (up to 
order 10) in the list of 500 frequencies. With the usual crite- 
rion that frequencies should be equal to within the square root 
of the sum of the squares of the individual 1-sigma errors, we 
find 983 such harmonics. However, simulations as in Sect. |2.5| 
show a factor 2-3 more harmonics. Our findings must therefore 
be attributed to coincidences, not to true harmonics. 

A similar search for linear combinations of the type f\+fi = 
/3 among the 500 frequencies reveals 45,412 results. But, again, 
simulations show these to be spurious. If we limit our search to 
frequencies significant under the red-noise hypothesis, we find 
no harmonics and only one linear combination. The reason we 
find so few of these is that most of our red-noise frequencies are 
between 10 and 20 d _1 ; harmonics and sums of these are then 
outside the 10 - 20 d _1 range. Because of the limited value of 
these results, we do not list them. 

In addition to these linear combinations, we investigate the 
frequency spacing for HD 46223. Degroote et al.| ( [20 10) found 
such spacing for the 08 V star HD 46149. Following their work, 
we look at both the autocorrelation function of the power pe- 
riodogram and the number of components we find in a chain 
(/(b fo + A/, /o + 2Af, ...) with an assumed frequency spacing 
A/. In both the chain and autocorrelation test, we start from the 
500 frequencies, but eliminate those outside the 3.0-12.0 d _1 
domain. The lower limit is taken to avoid the stronger peaks, 
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Fig. 6. Time-frequency analysis of the light curves. The top row uses a sliding window of 10 days, the bottom row one of 3 days. 
Left panels: HD 46223, middle: HD 46150, right: HD 46966. The grey scale indicates the semi-amplitude of the frequencies. It is 
clear that the frequencies do not have long lifetimes. 



the upper limit to avoid effects from the satellite orbital period 
(13.972 d~'). To determine the number of components in a chain 
at a given frequency, we take that /o + nAf equals a listed fre- 
quency if both differ by less than the step size we use in the 
figure, which is 0.015 d . For A/, we consider only the range 
of 0.1 to 2.0 d , as this is where both the autocorrelation and 
chain-test function have their largest values. 

Fig.B](top) shows the results for HD 46223. The autocorrela- 
tion function has a peak at A/ = 0.5 d _1 , but the number of com- 
ponents does not confirm this peak. We therefore conclude that 
we do not detect significant frequency spacings in HD 46223. 

2.11. Time-frequency analysis 

We also perform a time-frequency analysis (Fig. |6]l in order to 
detect variations in the amplitude of frequencies or in the fre- 
quencies themselves as a function of time. 

For this purpose, a sliding window with lengths of 10 and 
3 days is applied to the light curve and shifted with a step of 
one day. Then we apply a Fourier analysis on each frame by us- 
ing the HMM method (for further details on this time-frequency 
analysis method, see Huat et al. 2009 ). A plot of all frequencies 
in the 3-12 d _1 range detected on a sliding window of 10 and 
3 days is shown in Fig. [6] This figure confirms what we already 
found from the test where we split the data into two halves. The 
frequencies are clearly not stable over the ~ 34-day timescale of 
the observation. The figures do not show many frequencies that 
have a lifetime longer than the duration of the sliding window. 
This behaviour is compatible with the dominant presence of red 
noise. 

The limited lifetime also calls into question the applicabil- 
ity of the 1.5 IT criteri on (|Loumos & D eeming 1978) to resolve 
frequencies (see Sect. 2.5 i. A limited lifetime translates into a 



broadening of the frequency. This will reduce the number of 
unique frequencies found in our analysis. 



3. HD 46150 

For the analysis of HD 46150, we apply a similar detrend- 
ing as for HD 46223. The range of variations observed in the 



HD 46150 CoRoT light curve (Fig.|2j middle) is of order 16,000 
counts, corresponding to 8 mmag. Further details are listed in 
Table [T] No clear pattern in the variations is visible. The ampli- 
tude spectrum (Fig. [3] middle) is noisy and has its highest peak 
at 0.055 d _1 . Because there is little power at high frequencies, we 
limit the frequency domain to 0-100 d _1 in the further analysis. 

We next apply our standard prewhitening procedure and 
check the stopping criteria (Eq. (TV From the values found 
(Table[T|, we decide to limit the number of terms to 500. This list 
of frequencies and corresponding semi-amplitudes and phases is 
given in Table 3 (this table is only available in electronic form 
at CDS). Applying the significance criterion (Eq.|6]) with the as- 
sumption of white noise results in all 500 terms being signifi- 
cant. ThelLoumos & Deeming ( 1978 1 test shows 296 of these to 
be unique. Multifrequency fitting (see Sect. |2.6| l gives basically 
the same set of terms, though with slightly different amplitudes. 

However, it is again clear that there is a substantial red-noise 
component present in the spectrum, which we fit with the fit- 
function of Eq. [8] (see Table [T] and Fig. HI middle). Judging the 
significance of peaks with respect to this red noise, we find 50 
significant terms (indicated with 'R' in Table 3), using a 0.01 
cut-off. 

When we split the observed data set into two, we find that 
195 of the 500 terms occur in both halves, but this is less than 



the simulations of Sect. 2.9 show. Under the red-noise hypoth- 



esis, only seven of the 50 terms are present in both halves. We 
also explore the linear combinations of all 500 frequencies, but 
simulations show these to be not significant, as for HD 46223 
(see Sect. 2.10), Among the red-noise frequencies, we find no 



linear combinations, unless we relax the criterion we take for 
considering two frequencies to be equal. As for HD 46223, this 
shows that there are no significant linear combinations. 

We next apply a similar frequency spacing test as for 
HD 46223. A peak in the autocorrelation function is found at 
A/ = 0.58 d , but it is not confirmed by the plot with the num- 
ber of components in a chain (Fig. [5] middle). Finally, a time- 
frequency analysis (Fig. [6] middle) shows short "lifetimes" for 
the frequencies, which is again consistent with red noise. 
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Fig. 7. Evolutionary tracks relevant for our 
three stars. The models in the bottom figures in- 
clude overshooting, the ones at the top do not. 
Left: Tracks of 40 to 55 Mq computed with a 
mass-loss rate of 6.0 x 10~ 8 MQyr -1 . The er- 
ror boxes of HD 46223 (the hotter star) and 
HD 46150 (the cooler star) are also displayed. 
To show that both stars can be modelled by the 
same mass-loss rate without affecting the ex- 
cited frequency domain, we also plot a 49 Mq 
track computed with M = 6.8 x 10~ 8 MQyr~' 
(magenta-coloured line). Right: Same, but for 
HD 46966 computed with a mass-loss rate of 
5.0 x 10~ 9 MQyr _I . The full circles indicate po- 
sitions along the tracks that have excited fre- 
quencies. 



4. HD 46966 

The treatment of the data of HD 46966 is similar to that of 
HD 46223 and HD 46150 and details are listed in Table [Q The 
CoRoT curve of HD 46966 (Fig. [2] bottom) exhibits an almost 
regular oscillation with an approximate period of 10 days, but 
only about three cycles are covered by the observations. In ad- 
dition, we detect variability on shorter timescales. The peak-to- 
peak variation amplitudes correspond to about 6 mmag. 

Because there is little power at high frequencies, we limit 
the frequency domain to 0-100 d . We apply our standard 
prewhitening procedure and check the stopping criteria (Eq. T). 
From the values found (Table [TJ, we decide to limit the number 
of terms to 300. The list of frequencies and corresponding semi- 
amplitudes and phases is given in Table 4 (this table is only avail- 
able in electronic form at CDS). The |Loumos & D eeming ( 1978| l 
criterion shows 196 of these to be unique. The amplitude spec- 
trum (Fig. [3] bottom) is characterised by a first peak that is three 
times larger than the second one. This peak, at/ = 0.084 d _1 , i.e. 
12 days, corresponds to the observed almost regular oscillation. 

As for the other two stars, it is better to consider red noise, 
described by the fit-function from Eq. [8] Applying this, we 
find 53 significant terms, with a cut-off of 0.01. They are in- 
dicated with an 'R' in Table 4. Note that in this analysis, the 
/ = 0.084 d _1 frequency discussed above is no longer consid- 
ered significant. Table 4 shows that it has a false-alarm proba- 
bility of 0.997, which is too high compared to our 0.01 cut-off. 
Multifrequency fitting (see Sect. |2.6 1 gives basically the same 
set of frequencies, though with slightly different amplitudes. 

The test where we split the observed data set into two, shows 
that only 77 of the 300 terms are found in both halves. For the 
terms that are significant under the red-noise hypothesis, only 
three are found in both halves. Simulations as in Sect. |2.9| show 
that these numbers are not significant. 

We also explore the linear combinations between the 300 fre- 
quencies, but simulations show that these are not significant, as 
for the two other stars. Limiting ourselves to the 53 red-noise 
frequencies, we find no linear combinations. Only when we in- 



crease by a factor of three the criterion we take for equality of 
frequencies, do we start finding linear combinations. The result 
is therefore not significant. 

We also apply a similar frequency spacing test as for 
HD 46223. Again, a peak in the autocorrelation function is 
found, at A/ = 0.25 d _1 . The plot with the number of frequencies 
in a chain peaks at a different value of A/ = 0.43 d , however 
(Fig. [5] bottom), showing the result to be not significant. As for 
the two other stars, the time-frequency analysis (Fig. [6] bottom) 
shows short lifetimes for the frequencies; this again indicates the 
presence of red noise. 



5. Theoretical pulsations 

An important goal of the CoRoT photometry is to detect radial 
or non-radial pulsation modes in massive stars. To predict which 
modes can be expected, we compute a grid of models by using 
the Rome stellar evolution code ATON ( |Ventura et al.|2008] >. In 
these calculations, we assume a metallicity of Z — 0.015 and 
the metal mixture of Grevesse & Noels ( 1993 1. We also use the 
MAD non-adiabatic pulsation code ([Dupret et al. 2003 ) to de- 



termine the frequencies of excited modes. The computations are 
made on the basis of preliminary stellar and wind parameters 
derived from the atmosphere code CMFGEN (Hillier & Miller 
1998) by Martins et al. (in prep.). These parameters are listed 
in Table |5] We note however that we use the same value for the 
mass-loss rate for HD 46223 as for HD 46150 to compute the 
evolutionary tracks. The difference between the value used and 
those listed in Table [5] does not significantly affect the results 
mentioned in the present paper. Additionally, we consider both 
models with and without overshooting. 

HD 46223 and HD 46150 are in very similar positions in the 
HR diagram, whilst HD 46966 is cooler and more evolved. All 
three stars lie just at the boundary of the instability strip of the k- 
mechanism. Theoretical excited frequencies for that mechanism 
are generally in the range of 2-5 d" 1 for / = to I — 4 and a 
non-radial mode is also excited around 12 d~'. The theoretical 
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Fig. 8. Theoretical oscillation frequencies (for the /e-mechanism) 
computed with the MAD non-adiabatic code, as a function of 
the effective temperature. The results shown here do not in- 
clude overshooting. Top: 52 Mq track. Bottom: 30 Mq track. 
The black points indicate models that do not show pulsations, the 
coloured full circles stand for those models with the excited fre- 
quencies. The /-value is colour-coded (see legend). We also in- 
dicate, with vertical lines, the effective temperatures of the three 
stars studied here, as well as HD 46149 ( Degroote et al.| |2010) 
and HD 46202 ( [Briquet et al.|2011| >. 



Table 5. Stellar parameters for the three stars. 



Parameter 



HD 46223 



HD 46150 



HD 46966 



io g r eff (K) 


4.633 ± 0.020 


4.628 ± 0.021 


4.544 ± 0.025 


logg(cgs) 


4.00 


4.00 


3.75 


log L/Lq 


5.60 + 0.11 


5.65 + 0.11 


5.15 + 0.22 


R/R Q 


11.39 


12.35 


10.24 


logM(Moyr _1 ) 


-7.20 


-7.17 


-8.30 


v cq sin i (km s -1 ) 


100 


100 


50 


Vmacro (kms- 1 ) 


32 


37 


27 



Notes. Preliminary stellar and wind parameters derived by Martins et 
al. (in prep.). 



frequencies from our model calculations are shown on Fig.[8]for 
30 and 52 Mq (without overshooting): models with these masses 
present the largest number of excited modes. Models computed 
with overshooting give similar results. 



For HD 46223, we consider stellar evolution models with 
an initial mass ranging from 45 to 53 Mq (Fig. 17] left panel). 
On these tracks, models with a mass greater than 47 Mq show 
excited modes. For this mass range, the excited modes all have 
frequencies around 4 d _1 . However, if we look at the observed 
amplitude spectrum of HD 46223 (Fig. [3] top), no outstanding 
peak is situated close to 4 d~ ' which means that this star probably 
does not present pulsations, neither radial nor non-radial. For 
HD 46150, a similar conclusion can be drawn as for HD 46223. 
There are low-amplitude peaks around 4 d , but the peaks are 
not strong and they could be the result of noise. 

The amplitude spectrum of HD 46966 is clearly dominated 
by peaks at the frequencies / = 0.084 and 0.039 d -1 . The time 
coverage of the CoRoT observations is too short to know if these 
are regular variations. Furthermore, it is unlikely that these two 
frequencies are due to excited modes (Fig. IS] bottom) since those 
modes are predicted to be between 2.5 and 4 d -1 . Even though 
the amplitude spectrum has peaks situated in this frequency re- 
gion, these peaks are not strong. 

It is important to emphasize that the results presented in this 
section are dependent on the stellar parameters determined by 
Martins et al. (in prep.). These parameters are still under in- 
vestigation and should therefore be considered as preliminary. 
Nonetheless, the present work shows that even with the unprece- 
dented quality of the CoRoT light curves, we do not convinc- 
ingly detect the predicted excited modes. This suggests they are 
only weakly excited, or perhaps, not at all. On the other hand, a 
number of stronger peaks are observed, but at frequencies that 
are so different from the theoretical ones that we cannot at- 
tribute them to the expected excited modes, at least not within 
the framework of current theoretical models. It is thus not possi- 
ble to conclude with certainty that these three stars are pulsating. 
Moreover, we also explored the possibility of solar-like oscilla- 
tions, since Degroo te et al.| ((2010) found such oscillations in the 
OB binary system HD 46149. But the existence of such oscil- 
lations was rejected in HD 46150, HD 46223 and HD 46966 
because no regular spacings were found in the amplitude spec- 
tra. 



6. Discussion 

6.1. Radius and rotational frequency 

In a presumably single star, large-scale photometric variations 
can have a frequency compatible with the rotational frequency 
of the star. Their exact origin could be attributed to the existence 
of a spot on the star or be associated with the stellar wind. 

The first low frequency detected in the amplitude spectrum 
of HD 46223 is / = 0.249 d" 1 , corresponding to 4.02 days. 
With this period and the parameters from Table Bl we derive 
an equatorial velocity v eq of about 140 kms -1 which is greater 
than the v eq sin/ value of 100 kms _I 



measured by Mahy et al. 
(2009) and Martins et al. (in prep.). Consequently, the frequency 
/ = 0.249 d _1 could be associated with the rotational cycle of 
the star, implying a rotational period of about ~ 4 days and an 
inclination of 43 -45°. 

For HD 46150, the amplitude spectrum exhibits a strong 
peak at / = 0.055 d , corresponding to 18 days. The derived 
equatorial velocity is 34 kms -1 , which is too low compared to 
the v eq sin /. The frequency corresponding to the equatorial ve- 



locity found by |Mahy et al.| ( |2009] > or Martins et al., should be 
close to or larger than 0.16 d" 1 . The list of frequencies (Table 3) 
shows that the third frequency (/ = 0.144 d _1 ) in Table 3 is 
quite close. It corresponds to an equatorial velocity of 90 km s , 



10 



R. Blomme et al.: Variability in the CoRoT photometry of three hot O-type stars 



which is within acceptable range of the Ma hy et aL] and Martins 
et al. value, but would require an inclination of 90°. 

The amplitude spectrum of HD 46966 is marked by a strong 
peak at / = 0.084 d~\ representing a period of 12 days. With the 
stellar parameters mentioned in Table[5] we estimate the equato- 
rial velocity to be close to 44 km s . This velocity is just below 
the 1-sigma error bar on the v eq sin i determined by Mahy et al. 
(2009). Martins et al. (in prep.), on the other hand, clearly show 
that the v eq sin i is smaller than 50 km s . Consequently, the fre- 
quency / = 0.084 d _1 could be associated with the rotational 
cycle of the star, implying a rotational period of about 12 days. 



6.2. Pulsations and red noise 

For the three O-type stars studied here, a substantial number of 
sine function terms are found. However, it is unclear how many 
of those - if any - can be classified as true independent modes. 
They do not show the linear combinations or frequency spacings 
that may be expected in pulsators. If these terms were interpreted 
as pulsational modes, then their lifetimes would be very short. 
The amplitude spectra for these three stars are therefore unlike 
any of those of known pulsators in that part of the HR diagram 
( [Aerts et al.(20T0l Chapter 2). 

This is clearly illustrated by comparing the periodogram of 
these hot O stars (spectral types 04 - 08) to that of the cooler 
star HD 46202 (09 V). This star behaves like a classical pul- 
sator, with clear yS Cep-like pulsations that dominate the power 
spectrum (Briquet et al. 201 1 ). Various causes for what we see in 
our periodograms can be eliminated: rotational effects on pulsa- 
tion can be excluded as our periodograms are very different from 
those of Be stars (e.g. Diago et al. 2009). Spots on the surface 
of the star would give clearly isolated frequencies (and linear 
combinations between them). A power excess due to stochastic 
modes would occur at higher frequencies. 

The interpretation that ascribes most of the power in the pe- 
riodogram to red noise therefore seems more appropriate. We 
stress again that this red noise is caused by a physical mechanism 
in the star. It is not an instrumental effect, because its behaviour 
is different from star to star in these simultaneous observations. 
The parameters of the fit-function (TablefTli show that HD 46966 
has a somewhat lower a$ and a somewhat longer t than the other 
two stars. Compared with Plaskett's Star (HD 47129) a clearer 
difference is seen: especially the y value is different (y - 2.3, 
Mah y et "aI7||201 1 1. We also fitted the red-noise component of 
HP 46149 Pegroote et al."p0T0> . This star has a r = 29.5 ± 12 
and y - 0.55 ± 0.01, which are substantially different from the 
stars studied in this paper. All this points to an origin for the red 
noise that is intrinsic to the star. 

To see how convincing the red-noise description is, we sim- 
ulate light curves containing only red noise (according to the fit- 
function - Eq. [8]), and analyse them in the same way as the real 
data. We follow the prescription by Timmer & Koenig (1995) 
to generate the red noise. In the simulated data, the number of 
significant red-noise frequencies found is somewhat higher (by 
a factor less than two) than in the real data. Only one or two 
of these significant red-noise frequencies can be found in both 
halves of the light curve (somewhat less than for the real data). 
The number of linear combinations found is higher than for the 
real data, but this is presumably related to the higher number of 
significant red-noise frequencies. It therefore seems that the ob- 
served data can largely be described as red noise, though a few 
pulsational frequencies might be present. 



The presence of red noise in the stellar signal indicates 
stochastic, chaotic or quasiperiodic effects. It occurs in a num- 
ber of astrophysical contexts, such as the light curves of Mira 
variables (Templeton & Karovska 2009) and red supergiants 
(Kiss et al. 2006[), and t he X-ray light curves of active galax- 
ies ( Vaughan et al. 2003), dwarf novae ( Hakala et al.|2004| l and 
high-mass X-ray binaries ([Burderi 1994). 

The physical cause of this red noise in early-type stars is 
not clear. We next discuss three possibilities: sub-surface con- 
vection, granulation and inhomogeneities in the stellar wind. We 
stress that all three are highly speculative, and that as yet no 
modelling exists to see if they would give the correct quantita- 
tive behaviour. 

A link between red noise and convection was already pointed 
~ <[T975l> 



out by 



Schwarzschild 



Cantiello 



For early-type stars 
|et al.| ( |2009| ) show the existence of a sub-surface convection 
zone, caused by the opacity peak due to iron ionization. They 
show that stars in the HR diagram where this zone is present 
also have microturbulence, non-radial pulsations, wind clump- 
ing and line-profile variability, suggesting that these are caused 
by the convection zone. We speculate that the red noise detected 
in the early-type stars studied here is also a consequence of this 
convection zone, and could be associated with any of the above 
observational indicators. The position of our stars in the HR dia- 
gram shows that they should have a rather weak convection zone 
(see Cantiello et al. 2009 their Fig. 9). This suggests that the 
high quality of the CoRoT data allows us to see a more subtle 
effect than Cantiello et al. considered. Instead of attributing the 
red noise to convection only, there could also be an interaction 
between pulsation and convection. Perdang (2009) shows that 
this can result in red-noise dominated light curves. 

We also explore if the physical cause of the red noise could 



be described as granulation. Kallinger & Matthews d2010) show 
that, in the two 6 Scuti stars they study, the many hundreds of 
peaks in the amplitude spectrum can also be interpreted as being 
due to granulation. Their Fig. 3 relates the typical frequency of 
granulation (i.e. the inverse of the timescale) to a specific combi- 
nation of mass, radius and effective temperature. Although this 
relation is based on cooler stars, we nevertheless use it to pre- 
dict the granulation frequency in our stars. The value we find is 
~ 35 d , but this cannot be linked to any significant frequency 
or inverse timescale of the observations discussed in this paper. 
If granulation is the cause of the red noise in these hot O stars, 
it would therefore need to have properties that cannot simply be 
extrapolated from granulation in cooler type stars. 

Finally, the red noise could be due to inhomogeneities in the 
stellar wind. There are various observational indicators showing 
that the stellar winds of early-type stars are clumped (see re- 



view by Puis et al. 2008 , and references therein). The exact cause 



of this clumping is not clear. Intrinsic instabilities in the radia- 
tive driving mechanism have been proposed (Owocki & Rybicki 
1984), but the onset of clumping seems to happen very close to 
the stellar surface (Puis et al. 2006). The red noise could there- 
fore be related to the onset of clumping at the stellar surface. 

The behaviour of the early-type O stars (spectral types 04 
- 08) is considerably different from that of later-type stars. 
When taking out the dominant binary signature of HD 47129 
(08 III/I + 07.5 V/III), a series of frequencies is found that 
might be due to pulsations, but a clear red-noise component is 
present as well ( |Mahy et al."|20TT] ). HD 46149 (08, |Degroote 
|et al.|201~0| > does have significant frequency spacing. HD 46202 
(09 V) shows clear/? Cep-like pulsations (Briquet e t al.|20lT i. 
Observations with the MOST (Microvariability and Oscillations 
of Stars) satellite of the 09.5 V star f Oph also show /3 Cep-type 
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pulsations ( |Walker e t al. 2005). Amon g the early B-type stars, 
pulsation frequencies are easily found: Degroote et al. (2009a) 
present a study of 358 candidate B-star pulsators, showing nu- 
merous classical SPBs (Slowly-Pulsating B stars). We also note 
that a red-noise component is present in their spectra, but it is 
much more limited in frequency range than for the stars dis- 
cussed here (see their Fig. 2). 

The present data, combined with information from the liter- 
ature, permit a first attempt at empirically mapping the hottest 
part of the instability strip for /3 Cep and SPB stars. We can 
compare the position in the HR diagram of our stars, HD 47129 
dMahy et al |20TT), H P 46149 ( |Degroote et al.|2010j i, HD 46202 
( |Briquetetal.|2011| ) and t, Oph \ Walker et al.|200"5| t o theoreti- 
cal predictions of the /3 Cep and SPB instability strip. Zdravkov 
& Pamyatnykh (2008 ) present the theoretical strip for models up 



to 40 Mq for a metallicity Z = 0.02. Most of the stars in our 
list fall into their theoretical /3 Cep instability strip, but outside 
the SPB one. As only a few of these stars actually show f3 Cep 
pulsations, this suggests either that the strip is much narrower, 
or that the pulsations are excited for only a small fraction of the 
stars in the strip. Especially constraining is a group of three stars 
(HD 46966, HD 46202 and ( Oph) which are close to one an- 
other in the HR diagram, with HD 46202 and ( Oph showing /3 
Cep pulsations, but HD 46966 not. 

A number of O-type /J Cep pulsators have also been detected 
in ground-based observations (e.g. |Telting et al. 2006; D e Cat| 
et al. 2007; Pigulskf&Tojmanski 2008]). These are all late O- 
type stars and therefore close to the three-star group. A detailed 
knowledge of their astrophysical parameters would allow a bet- 
ter mapping of the instability strip in that part of the HR diagram. 
In the hotter part, we note the absence of/3 Cep pulsations in the 
main-sequence stars discussed here. The theoretical /3 Cep strip 
extends to the terminal age main sequence (TAMS) however, and 
stars more evolved than those studied here could therefore show 
P Cep pulsations. 

7. Conclusions 

The CoRoT light curves of the hot O-type stars HD 46223, 
HD 46150 and HD 46966 were analysed using standard meth- 
ods to search for pulsation frequencies. The detection of such 
frequencies would allow an asteroseismological interpretation. 
However, the results show that most of the variations of these 
stars are of a stochastic nature. The only significant exception to 
the above is the possible rotation period detected for HD 46223 
and HD 46966, although these conclusions are not strong and 
still require confirmation. Evidence for a rotation period in the 
HD 46150 data is even less strong. 

The periodogram of these three stars is clearly different from 
that of a classical pulsator such as HD 46202. Most of the power 
in our periodograms is not due to pulsation, but is more correctly 
described by red noise. Based on the limited number of stars 
studied so far, a trend is suggested: the earliest O stars show red 
noise, while the later O-types have pulsational frequencies of the 
ft Cep type. The switch-over occurs around spectral type 08. 

The specific physical cause of this noise in O-type stars is 
at present unclear. We point out the possibilities of sub-surface 
convection, granulation or inhomogeneities in the wind. All 
three options are highly speculative and therefore await confir- 
mation by detailed modelling. The CoRoT data of these three 
hot stars present interesting challenges for our understanding of 
the outer photosphere and inner wind of early-type stars. 
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